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INTRODUCTION 


J. 


Achieving  high  specific  impulses  by  increasing  the  solids  loading 
of  highly  energetic  materials  (e.g.,  HMX)  in  solid  propellants  has  been 
an  objective  of  the  rocket  industry  for  many  years.  At  a very  high 
solids  loading  there  is  an  increased  probability  that  cracks  or  flaws 
can  develop  in  the  propellant  grain.  These  defects  may  originate  in 
solid  propellants  as  a result  of  (1)  improper  manufacturing  processes  so 
that  bubbles  and  cavities  are  present  in  the  propellant,  (2)  the  aging 
effect  during  storage,  (3)  nonuniform  stress  loading  associated  with 
large  pressure  gradient  during  propellant  ignition  process,  (4)  the 
impairment  of  the  adhesive  bond  between  the  fuel  and  oxidizer  due  to 
mechanical  stresses  during  combustion,  (5)  the  formation  of  micropores 
due  to  thermal  stresses,  and  (6)  the  damage  of  propellant  charges  during 
transportation. 

On  a microscopic  level  it  is  believed  that  cracks  occur  in  hetero- 
geneous or  composite  propellants  due  to  an  adhesive  failure  between  the 
rubbery  binder  and  the  high  solids  loading  (up  to  approximately  90%  on  a 
volume  basis)  of  oxidizer  particles  and  other  granular  fuel  additives. 
Once  the  adhesion  between  these  particles  has  failed,  microscopic 
vacuoles  form  near  the  particles  and  readily  propagate  to  become  finite 
cracks.  In  general,  these  finite  cracks  can  degrade  the  mechanical 
properties  and  structural  integrity  of  solid  propellant  grains. 

The  extent  to  which  crack  combustion  can  reduce  the  reliability  of 
solid  rocket  motors  has  been  one  of  the  main  concerns  in  the  development 
of  high-energy  propellants.  The  cracks  can  allow  the  hot  gases  to 
penetrate  the  cavity  and  can  provide  additional  surface  area  for  ignition 
and  combustion.  The  combustion  processes  inside  a propellant  crack  can 
produce  shock-like  compression  waves  and  may  also  lead  to  a pressure 
build-up  which  could  cause  mechanical  deformation,  crack  propagation,  and 
possible  def lagration-to-detonation  transition  (DDT).  Even  if,  under 
certain  conditions,  the  convective  burning  does  not  lead  to  catastrophic 
failure,  the  combustion  inside  the  propellant  crack  can  cause  the  per- 
formance of  a rocket  motor  to  deviate  significantly  from  expectations. 

Although  imperfections  in  solid  propellant  grains  may  be  distributed 
in  an  irregular  intricate  network,  a rational  starting-point  for  under- 
standing the  phenomena  involved  in  crack  combustion  is  to  consider  a 
simple  ordered  system — a single  isolated  crack.  For  such  a system  it  is 
desirable  to  establish  a criterion  for  stable  combustion  without  cata- 
strophic crack  propagation  or  incipient  detonation;  the  delineation  of  a 
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safe  domain  of  governing  parameters  Is  the  ultimate  goal.  But  before  it 
can  be  achieved,  a realistic  theoretical  model  describing  both  the 
mechanical  deformation  and  transient  combustion  processes  inside  a 
propellant  crack  has  to  be  established  and  experimentally  verified. 

The  primary  objectives  of  the  research  effort  were  as  follows: 

1.  To  conduct  experimental  firings  with  single-pore  propellant 
grains  for  verification  of  the  theoretical  model  and  for  establishment 
of  a fundamental  data  base  in  crack  combustion. 

2.  To  modify  the  theoretical  model  and  computer  program  for 
obtaining  the  in-depth  thermal  profile  of  the  propellant  so  that  the 
dynamic  burning  behavior  can  be  studied. 

3.  To  consider  various  ignition  criterions  such  as  the  "Dipstick" 
criterion;  the  critical  surface  temperature  criterion;  or  the  two- 
temperature,  surface- ignition  criterion  for  the  calculation  of  flame 
spreading  rates. 

4.  To  nondimensionalize  the  gas-dynamic  and  heat-transfer  equa- 
tions in  order  to  obtain  a set  of  dimensionless  parameters  that  charac- 
terize crack  combustion. 

5.  To  study  the  compatibility  of  the  crack  combustion  program 
with  a nonlinear  stress  analysis  program  so  that  the  stress  field  near 
the  crack  tip  can  be  evaluated  for  determination  of  crack  propagation 
due  to  pressure  loading.  Also,  to  select  a well-developed  structural 
analysis  program  suitable  for  coupling  with  the  crack  combustion  code. 


BACKGROUND 


During  previous  research  investigations  conducted  at  the 
Pennsylvania  State  University,  substantial  progress  was  achieved  in  the 
study  of  the  convective  burning  processes  inside  a solid  propellant 
crack.1-3  Some  of  the  important  theoretical  and  experimental  achieve- 
ments are  given  as  follows: 


1 Naval  Weapons  Center.  Transient  Combustion  in  Solid  Propellant 
Cracks,  by  K.  K.  Kuo,  et  al. , Pennsylvania  State  University.  China 
Lake,  CA,  NWC,  October  1977.  (NWC  TP  5943,  publication  UNCLASSIFIED.) 

ZK.  K.  Kuo,  A.  T.  Chen,  and  T.  R.  Davis.  "Transient  Flame  Spread- 
ing and  Combustion  Processes  Inside  a Solid  Propellant  Crack,"  presented 
at  AIAA  15th  Aerospace  Sciences  Meeting,  Los  Angeles,  CA,  January  1977. 
(AIAA  Preprint  No.  77-14.) 

3K.  K.  Kuo  and  D.  R.  McClure.  "Transient  Combustion  Processes  in 
Solid  Propellant  Cracks,"  Semi-Annual  Report  to  the  Aerothermochemistry 
Division  of  the  Naval  Weapons  Center,  August  1977.  (Report  UNCLASSIFIED.) 
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1.  An  extensive  literature  survey  of  theoretical  and  experimental 
investigations  on  flame  propagation  in  cracks  has  been  conducted.  A 
survey  of  heat-transfer  and  flow-resistance  empirical  correJations  on 
flow  through  chemically  inert  cracks  was  performed. 

2.  A quasi-one-dimensional  theoretical  model  has  been  developed 
to  describe  the  transient  combustion  phenomena  in  a propellant  crack. 

The  theoretical  model  can  be  used  to  predict  the  wave  phenomena,  heat- 
transfer  processes  from  gas  to  propellant  surface  flame  propagation,  and 
resultant  pressurization  at  various  locations  along  the  propellant 
cavity . 

3.  A computer  program  has  been  developed  based  on  the  above 
theoretical  model  to  predict  the  flame-spreading  and  combustion  phenomena 
inside  a crack.  The  program  can  be  utilized  to  conduct  parametric 
studies  on  the  effects  of  crack  gap  width,  depth,  cross-sectional  area 
distribution,  propellant  burning  rate,  ignition  (runaway)  temperature, 
and  other  phys iochemical  properties  of  the  propellant. 

4.  A laboratory-sized  solid-propellant  combustion  chamber,  util- 
izing current  state-of-the-art  techniques  in  high-speed  photography  and 
transient  pressure  recording,  has  been  developed  to  obtain  data  for 
verification  of  the  theoretical  model.  The  facility  permits  the  simul- 
taneous acquisition  of  the  time-correlated  flame-spreading  rate  and 
pressure  wave  propagation. 


THEORETICAL  WORK 


NONDIMENSIONAL IZATION  OF  THE  GAS -DYNAMIC  AND  HEAT- TRANSFER  EQUATIONS 
GOVERNING  THE  CONVECTIVE  BURNING  PROCESSES  IN  CRACKS 

The  Russian  work4*4  suggests  that  the  transition  from  conductive 
(normal  layer-to-layer)  burning  to  convective  (gas-penetration)  burning 
for  gas-permeable  propellants  occurs  when  the  Andreev  number.  An,  is 
above  a critical  value.  The  Andreev  number  is  defined  as 


4A.  F.  Belyaev,  et  al.  "Development  of  Burning  in  a Single  Pore," 
in  Transition  of  Combustion  of  Condensed  Systems  to  Detonation . Science 
Publisher,  Moscow,  1973.  Chapter  3,  Part  B and  Chapter  5,  Part  A. 

4A.  D.  Margolin  and  S.  V,  Chuiko.  "Combustion  Instability  of  a 
Porous  Charge  with  Spontaneous  Penetration  of  the  Combustion  Products 
into  the  Pores,"  Combustion , Explosion  and  Shock  Waves.  Vol.  2,  No.  3, 
1966. 
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An  = 


p r C d 
c b pr  p 


X 

Pr 


(1) 


where 


p = density  of  the  gas-permeable  (porous)  charge 
r^  = burning  rate 

Cpr  = specific  heat  of  the  condensed  phase 
d = diameter  of  the  pore 


P _ 


= thermal  conductivity  of  the  propellant 


The  Andreev  number  is  similar  to  the  Peclet  number  in  its  physical 
meaning;  it  represents  the  ratio  of  convective  heat  flux  to  conductive 
heat  flux.  This  finding  is  helpful  in  defining  the  boundary  beyond 
which  there  is  spontaneous  penetration  of  the  combustion  products  into 
the  pores. 


Once  the  gas  has  penetrated  into  the  pores  of  a gas-permeable 
propellant,  the  Andreev  number  no  longer  governs  the  flame  spreading  and 
combustion  phenomena  inside  the  pore.  Several  questions  can  be  raised 
immediately  in  regard  to  the  governing  dimensionless  parameters  which 
dominate  the  physical  processes  in  propellant  crack  combustion.  First, 
what  are  the  dimensionless  parameters  which  are  important  to  consider  in 
the  convective  burning  stage  of  the  multifaceted  problem  of  DDT?  Second, 
is  the  Andreev  number  the  only  parameter  in  governing  the  gas  pene- 
tration condition  into  cracks?  Would  the  rate  of  gas  pressurization  at 
the  crack  opening  also  be  important  in  determining  the  state  of  gas 
penetration  at  the  crack  entrance?  Third,  would  the  Andreev  number  be 
a natural  dimensionless  group  obtained  from  the  governing  equations?  If 
not,  can  it  be  constructed  from  other  dimensionless  parameters  obtained 
from  the  nondimensionalization  procedure? 


To  answer  these  questions  and  also  to  determine  the  controlling 
dimensionless  parameters  that  characterize  the  convective  burning  prob- 
lem, it  is  worthwhile  to  nondimensionalize  the  governing  gas-dynamic  and 
heat  transfer  equations.  It  should  be  noted  that  the  complete  groups  of 
governing  parameters  cannot  be  obtained  without  looking  at  the  coupled 
solid  mechanics  equations  for  the  solid  propellant.  However,  for  the 
time  being,  it  is  useful  to  obtain  a set  of  dimensionless  parameters 
from  the  crack  combustion  model  before  coupling  it  to  solid  mechanics; 
the  controlling  parameters  may  suggest  a systematic  way  of  performing 
laboratory  tests  and  logical  ways  of  correlating  experimental  data. 
Therefore,  the  following  nondimensionalization  procedures  are  carried 
out. 

To  simplify  the  nondimensionalization  procedure,  the  terms  of 
secondary  importance  in  the  gas-dynamic  equations  were  neglected.  After 
neglecting  the  body  force  term  and  the  normal  viscous  stress  gradient 
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Some  of  the  reference  parameters,  designated  with  superscript  *,  can  be 
selected  readily  without  any  difficulty.  For  example,  the  character- 
istic length  along  the  crack  should  be  the  length  of  the  crack,  L,  i.e.. 


* 

x = L 


The  reference  port  area  can  be  taken  to  be  the  initial  average  port  area 

°H  Pw 
. _ i i 

\ = 4 (15) 

where  Dh^  represents  the  initial  hydraulic  diameter  of  the  pore.  The 
reference  burning  area  can  be  selected  to  be 

\ " 4'b1  (16) 


H * *■ 
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Naturally,  the  initial  average  wetted  perimeter  pw^  and  the  initial  aver- 
age burning  perimeter  p should  be  used  a.  reference  quantities  shown 
in  Eq.  (11).  bl 

Based  upon  the  definition  given  in  Eq.  (8),  the  reference  temper- 
ature for  the  gas  temperature  is  selected  to  be 

T*  Tf  - Tt  (17) 

The  reference  gas  density  can  be  defined  by  the  perfect  gas  law  as 


•k  p 

P - — (18) 

RT 

* 

where  P , the  reference  pressure,  must  be  properly  selected;  and  the 

value  of  P*  can  be  used  to  determine  the  burning  rate,  r * (.a 

b 

rb  = a(P  )n  (19) 

After  balancing  the  mass  burning  rate  with  the  mass  flow  rate  out  of  a 
crack  at  the  choked,  condition,  P*  can  be  selected  as 


(20) 


(21) 


The  most  difficult  reference  parameter  to  select  in  the  nondimen- 
sionalization  is  the  appropriate  t*.  There  are  many  characteristic 
times  one  can  consider  in  the  crack  combustion  study.  To  name  a few,  we 
have:  flow  residence  time,  flame  propagation  time,  ignition  delay  time, 

thermal  wave  relaxation  time  in  solid  propellant  (time  associated  with 
thermal  wave  propagation  during  propellant  regression),  pressure  wave 
traveling  time  from  crack  entrance  to  tip,  characteristic  time  related 
to  the  initial  pressurization  rate  at  the  crack  entrance,  the  time  to 
reach  maximum  pressure  in  the  overall  transient  event,  etc.  These  times 
are  not  of  the  same  order  of  magnitude.  Also,  these  reference  times  are 
associated  with  different  length  scales  and  velocities.  There  are  many 
length  scales  that  can  be  enumerated.  For  example,  crack  length,  gap 
width,  hydraulic  diameter,  crack  surface  roughness,  flame  stand-off 
distance,  thermal  wave  penetration  depth  due  to  convective  heating, 
heated  layer  thickness  of  the  burning  propellant,  etc. 
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After  a detailed  consideration,  the  ignition  delay  time,  which  is 
closely  related  to  the  flame  spreading  time,  is  taken  to  be  the  most 
important  time  scale  in  convective  burning.  By  setting  Tp^- = Tj  and 
nondimensionalizing  Eq.  (5)  according  to  the  definitions  given  in 
Eq.  (6),  (8),  (9),  and  (13),  we  have 


* 

where  h , the  reference  heat-transfer  coefficient,  must  be  defined. 
Before  £he  explicit  form  of  h * is  defined,  it  is  interesting  to  see 
the  physical  meaning  of  the  following  parameter  formed  from  Eq . (23) 


* 

t a h 
pr  c 


Pr 


* 2 

(Bi  r 


(24) 


where  Bi  , the  Biot  modulus,  represents  the  ratio  of  internal  thermal, 
resistance  of  the  propellant  to  the  external  thermal  resistance  of  the 
convective  boundary  at  the  gas-propellant  interface. 


Based  upon  the  Reynolds  analogy  for  turbulent  flow  in  a tube  and 
the  well-established  Blasius  flow-resistance  formula,  h * can  be  related 

it  C 

to  u according  to 

h * = 0.0038  C pV  (25) 

c p 

it  it 

The  characteristic  velocity,  u , can  be  related  to  t by 


(26) 


The  advantage  of  using  the  above  relationship  between  u and  t is  to 
have  both  the  time-variation  and  the  space-variation  of  the  mass,  momentum 
and  energy  in  the  governing  equations  in  the  same  order  or  magnitude. 


10 


NWC  TP  6049 


In  addition  to  the  preceding  reference  parameters,  we  further 
define 

M*  * 8.699  x 10~7  x (Mw)°‘5  (T*)0-65  ' (27) 

based  upon  Svehla*s  formula.6  The  reference  wall  shear  stress  on  crack 
surfaces  can  be  defined  by  using  the  Weisbach  friction  coefficient,  f, 
as 


* 

T 

W 


f * *2 

8g  P U 


Also,  we  define 


* is 

E = C T and 
P 


* 

M 


*2 

u 


YgRT* 


(28) 


(29) 


After  nondimensionalization,  the  governing  gas-dynamic  equations 
(2)-(4)  become 


6National  Aeronautics  and  Space  Administration.  Estimated  Vis- 
cosities and  Thermal  Conductivities  of  Gases  at  High  Temperatures , 
by  R.  A.  Svehla.  Lewis  Research  Center,  Cleveland,  OH,  NASA,  1962. 
(NASA  TR  R-132,  publication  UNCLASSIFIED.) 
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The  dimensionless  parameters  shown  in  the  governing  equations  (30)-(32) 
and  (22)  are  [through  the  use  of  Eq.  (23)  and  (24)]: 


* * ’ A • A 

P u PA  P± 


f, 


* 

M 


* 

h 

c 

* * 

p u C 

P 


* 

Bi 


The  first  group  can  be  defined  as  the  mass  blowing  parameter,  B 

m 

* 

g * ^P^k  _ mass  burning  rate  per  unit  burning  area 
m * * mass  flow  rate  per  unit  port  area 


* * 

P u 


(33 


The  second  group  can  be  expressed  as 


(34) 


When  the  burning  perimeter  is  equal  to  the  wetted  pewmeter,  the  area 
ratio  is  directly  related  to  the  length-to-hydraulic  diameter  ratio  of 
the  crack.  The  third  group  can  be  replaced  by  4(L/D  ). 

Hi 

The  friction  coefficient  for  turbulent  flow  Inside  a crack  of  a 
given  surface  roughness  and  condition  of  gasification  can  be  expressed 
in  the  following  functional  form 
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f *=  f (Re 


* 

Bm  > 

m 


(35) 


where  the  Reynolds  number.  Re  , can  be  defined  by 

* * 

* P u °H 

Re  D £ i 

H.  * 

i P 


(36) 


and  e/Dn^  represents  the  ratio  of  surface  roughness  to  the  hydraulic 
diameter  of  the  crack.  The  ratio  of  h * to  p*u*C  can  be  defined  as 


the  Stanton  number,  St  , i.e. 

* 


St 


* * 

p u C 


(37) 


Other  dimensionless  parameters  which  should  be  considered  are: 


Among  these  parameters,  the  first  group  comes  from  the  boundary  condi- 
tions at  the  crack  entrance;  it  represents  the  ratio  of  the  thermal  wave 
relaxation  time  for  burning  propellant  to  the  time  characterizing  the 
initial  pressurization  at  the  crack  entrance.  The  second  group  represents 
the  ratio  of  the  flame  stand-off  distance  to  the  crack  gap  width;  it  is 
important  in  the  consideration  of  incomplete  combustion  due  to  an  extremely 
small  gap  width  between  the  two  opposite  crack  surfaces.  The  last  group 
represents  the  ratio  of  the  unknown  flame-propagation  rate  to  the  char- 
acteristic velocity  u*  defined  in  Eq.  (26). 


Finally,  these  dimensionless  parameters  can  be  put  into  the  follow- 
ing functional  form 


fs 


gap 


ft  n 

Bi  » B„ 


* * 

M , St  , 


T.  -T. 
jgB-J: 


VTi 


(38) 
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It  should  be  noted  that  these  parameters  do  not  include  the  effect  of 
mechanical  deformation  and  damage.  In  our  continuing  research  effort  on 
crack  combustion,  the  governing  parameters  associated  with  solid  mechan- 
ics will  be  studied  and  incorporated  into  the  dimensionless'  group . 


SOLUTION  OF  THE  SUBSURFACE  TEMPERATURE  DISTRIBUTION  FROM  TRANSIENT 
HEAT-CONDUCTION  EQUATION 


There  are  two  major  reasons  for  solving  the  transient  heat- 
conduction  equation.  One  reason  is  to  incorporate  the  ignition  criter- 
ion, such  as  the  "Dipstick"  ignition  criterion,  which  requires  the 
information  of  subsurface  temperature  profile.  The  other  reason  is  to 
facilitate  future  investigation  of  the  dynamic  burning  effect  in  propel- 
lant crack  combustion.  Recall  that  previously  the  propellant  surface 
temperature  at  various  stations  along  the  crack  wall  was  determined  by 
Goodman's  integral  method7  described  in  footnote  1.  Currently,  in  order 
to  obtain  the  information  about  the  in-depth  thermal  response  of  the 
propellant,  the  transient  heat  conduction  equation 


with  the 


3T  (y , t)  32T  (y , t) 

pr  J m pr  7 

3t  upr  2 

9y 

following  initial  and  boundary  conditions 


(39) 


T 

pr 


(y.0) 


(40) 


V-*) 


(41) 


3Tpr(°,t) 

3y 


hc(t) 

X 


pr 


[T 

gas 


(t) 


Tpr(°,t)] 


(42) 


has  been  coded  in  a new  subroutine  for  solution  by  an  explicit,  variable- 
mesh,  finite-difference  method.  Equation  (39)  is  used  to  obtain  the 
propellant  temperature  distributions  prior  to  ignition.  The  finite- 
difference  and  Goodman  integral  methods  have  been  compared  to  an  exact 
solution8  of  the  heat  flow  in  a semi-infinite  solid  with  constant  h . 

Both  approximate  methods  were  found  to  be  quite  accurate  (<3%)  . The  in- 
depth  temperature  profiles  will  be  used  in  future  calculations  to  eval- 
uate the  effects  of  dynamic  burning  and  the  importance  of  thermal  stresses 
in  the  finite  element  work. 


7T.  R.  Goodman.  "The  Heating  of  Slabs  with  Arbitrary  Heat  Input," 
Journal  of  the  Aerospace  Sciences,  Vol.  26,  March  1959,  pp.  183-188. 

8P.  J.  Schneider.  Conduction  Heat  Transfer.  Addison-Wesley 
Publishing  Company,  Reading,  Mass.,  1955. 
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The  "Dipstick"  ignition  criterion,  postulated  by  Price,  has  been 
incorporated  in  the  crack  combustion  program.  The  concept  of  a critical 
gasification  temperature  and  ignition  temperature  was  adde<l.  to  this 
criterion  to  provide  easier  interfacing  with  the  previous  two- 
temperature,  surface-ignition  criterion  (Tprcr^  — Tprign  = 460°K  at 

^ign  * 60  pm)  which  revealed  that  the  propellant  surface  temperature 

near  the  crack  entrance  may  exceed  1000-1300°K  before  the  substrate 

temperature  at  6jgn  reaches  the  ignition  temperature.  However,  the 

onset  of  ignition  can  be  established  earlier  by  lowering  T or 

reducing  6.  . prcri 

ign 

To  have  the  maximum  flexibility,  the  previously  developed  sub- 
routine which  uses  the  integral  solution  for  the  propellant  surface 
temperature  calculation,  is  retained  in  the  computer  program  as  an 
optional  method  in  the  computation  of  the  flame  spreading  rate. 


1 


IMPROVEMENTS  IN  THE  NUMERICAL  TREATMENT  OF 
GAS-PHASE  BOUNDARY  CONDITION 

* 

As  pointed  out  by  Price,  some  steps  in  the  numerical  treatment  of 
the  gas-dynamic  boundary  condition  for  inflow  could  be  improved.  This 
deficiency  has  been  rectified,  and  the  program  can  now  more  accurately 
describe  the  influence  of  rapidly  varying  rocket  chamber  pressurizations 
on  the  subsonic  and  sonic  inflow  of  gases.  Some  details  of  these 
changes  are  discussed  below. 

In  the  numerical  computation  of  the  gas  dynamic  boundary  condition 
for  the  case  of  inflow  into  the  crack,  it  was  found  that  a quasi-steady 
approximation  can  be  used  to  treat  the  flow  from  the  rocket  chamber  into 
the  crack  entrance  region.  As  long  as  the  boundary  station  for  partial 
differential  equation  solution  is  taken  to  be  very  close  to  the  actual 
crack  opening  station  (AL  £ 0.01  cm),  the  pressure  and  gas  temperature 
at  the  boundary  station  can  be  accurately  determined  from  the  isentropic 
relations  of  the  actual  or  simulated  rocket  chamber  conditions.  The  gas 
velocity  at  the  boundary  station  is  determined  from  an  integration  of 
the  particle  pathline  in  the  Left  Boundary  Condition  Subroutine  (LBC) . 
When  the  pressure  at  the  boundary  station  is  greater  than  the  critical 
pressure  required  for  choked  flow,  the  gas  velocity  from  the  previous 
iteration  and  the  new  values  of  external  stagnation  pressure  and  temper- 
ature are  used  in  a subsonic,  isentropic  relation  to  evaluate  the  new 
pressure  and  temperature  at  the  boundary  station.  It  is  found  that  this 
technique  works  very  well  even  for  a rapidly  varying  rocket  chamber 
pressure . 


Private  communication  with  Channon  F.  Price,  Naval  Weapons  Center, 
April  1977. 
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STUDY  OF  STRUCTURAL  ANALYSIS  PROGRAM 

A literature  survey  on  crack  fracture  and  deformation  has  been 
conducted.  The  review  of  various  well-developed  finite  element  analysis 
programs  has  also  been  carried  out.  Of  the  numerous  structural  analysis 
codes  (e.g.,  NONSAP,  NASTRAN,  TEXGAP) , the  modified  version  of  NONSAP 
(NFAP)  is  considered  to  be  the  most  suitable  one  to  couple  with  the 
crack  combustion  code  for  the  reasons  discussed  below. g_1 1 

The  modified  version  of  NONSAP  (NFAP)  has  been  obtained  from  the 
University  of  Akron,  Akron,  Ohio  and  loaded  on  the  PSU  computer.  The 
user's  manual  has  been  reviewed  to  determine  the  best  way  to  streamline 
and  couple  the  program  with  the  combustion  code.  The  important  features 
of  NFAP  which  are  relevant  to  the  research  work  on  crack  combustion  are: 

1.  The  program  is  an  efficient,  easily  modifiable  general  stress 
analysis  program  with  restart  capabilities. 

2.  It  can  treat  transient  loading  and  inertia  effects  in  two-  and 
three-dimensional  structures. 

3.  The  code  can  treat  geometric  nonlinearities  caused  by  large 
deformations. 

4.  NFAP  also  can  treat  material  nonlinearities.  It  contains  a 
special  dummy  overlay  in  which  the  user  can  supply,  in  a consistent 
manner,  one's  own  form  of  the  propellant  constitutive  equation. 


EXPERIMENTAL  SET-UP  OF  THE  TEST  RIG 


CRACK  COMBUSTION  CHAMBER 

The  major  components  of  the  crack  combustion  chamber  used  during 
this  reporting  period  are  essentially  the  same  as  those  described  pre- 
viously in  the  last  annual  report.1  The  exit  nozzle,  which  controls  the 
peak  chamber  pressure  and  affects  the  chamber  pressurization  rate,  were 

9K.  J.  Bathe,  E.  L.  Wilson,  and  R.  H.  Iding.  NONSAP:  A Struc- 
tural Analysis  Program  for  Static  and  Dynamic  Response  of  Nonlinear 
Systems.  University  of  California,  Berkely,  CA,  February  1974. 

(USESM  74-3,  publication  UNCLASSIFIED.)  (Also  NTIS  PB-231  112.) 

10T.  Y.  Chang,  S.  Prachuktam,  and  M.  Reich.  "Assessment  of  a Non- 
linear Structural  Analysis  Finite  Element  Program  (NONSAP)  for  Inelastic 
Analysis,"  presented  at  ASME-Energy  Technology  Conference,  18-23  Septem- 
ber 1977,  Houston,  Texas.  (ASME  Paper  77-PVP-10.) 

1 'T.  Y.  Chang  and  S.  Prachuktam.  NFAP,  A Nonlinear  Finite- Element 
Analysis  Program.  University  of  Akron,  Akron,  Ohio,  October  1976. 
(Report  No.  SE  76-3,  publication  UNCLACSIFIED. ) 
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previously  made  of  stainless  steel.  After  several  test  firings,  the 

exit  nozzle  throat  area  was  enlarged  due  to  metal  erosion;  usually,  the 

area  became  noncircular,  which  introduces  difficulties  in  achieving 

reproducible  test  runs.  The  metal  erosion  problem  has  been  alleviated 

by  fabricating  the  nozzle  from  titanium  which  is  less  sensitive  to  metal  l 

erosion.  The  titanium  nozzle  has  shown  no  metal  erosion  in  test  firings 

done  to  date. 


IGNITION  SYSTEM 

There  are  two  major  drawbacks  of  the  previously  used  methane-oxygen 
igniter.  First,  the  gaseous  ignition  system  produces  ignition  with  a 
long  ignition  delay  time  followed  by  a very  short  flame  spreading  inter- 
val. In  order  to  observe  the  flame  spreading  phenomena  along  the  crack 
with  a proper  framing  rate,  one  must  either  use  a very  high  footage  of  ' 

film  or  be  able  to  control  the  Initial  starting  of  the  camera  to  incor- 
porate the  long  ignition  delay.  The  use  of  high  footage  films  is  ex- 
tremely expensive  and  the  control  of  the  camera  start  time  would  be  very 
difficult  to  achieve.  Second,  as  described  in  Ref.  1,  the  rocket  chamber 
pressurization  rate  was  found  to  be  an  important  variable  in  influencing 
the  flame  propagation  rate  inside  propellant  cracks.  The  gaseous  igniter 
can  produce  only  relatively  low  chamber  pressurization  rates.  In  order 
to  simulate  the  conventional  pressuriztion  rates  during  rocket  motor 
ignition  transient  or  to  simulate  the  case  of  submerged  cracks,  it  is 
desirable  for  the  igniter  to  be  capable  of  obtaining  broad  ranges  of 
chamber  pressurization  rates. 

Based  on  these  requirements,  it  was  decided  to  develop  a more 
versatile  ignition  system  than  the  previously  used  gaseous  ignition 
system.  Prior  to  the  development  of  a successful  ignition  system,  a 
variety  of  igniter  designs  were  tested  and  found  to  be  inadequate  during 
the  early  development  phase. 

Numerous  experiments  were  conducted  with  an  igniter  which  used  a Mk 
125  ignition  element.  This  igniter  produced  an  undesirable  amount  of 
metallic  debris  and  generated  uncontrollably  large  pressure  excursions 
in  the  small  volume  near  the  entrance  to  the  propellant  crack.  Conse- 
quently, the  Mk  125  igniter  design  was  abandoned,  and  a different  igni- 
tion system  was  designed  which  would  permit  testing  with  a variety  of 
propellant  charges. 

As  schematically  described  in  Figure  1,  the  igniter  uses  a solenoid- 
actuated  firing  pin  to  initiate  a percussion  element  (FA-34)  which 
ignites  propellant  shavings  in  the  booster  chamber  and  the  propellant  in 
the  main  igniter  chamber.  A variety  of  gun  propellants  (WC  870,  WC  846, 

M 6,  and  M 26)  and  solid  rocket  propellants  (AFRPL  nonmetallized  AP 
composite  and  Navy  Propellant  Type  A)  were  tested  in  the  main  igniter 
chambers.  Only  the  rocket  propellants  produced  a workable  system  that 
could  be  easily  ignited  and  would  not  prematurely  extinguish. 
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FIGURE  1.  Schematic  Diagram  of  Propellant  Crack 
Combuation  with  the  Solid  Propellant  Igniter  System. 
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Although  the  propellant  shavings  were  useful  in  helping  the  igni- 
tion of  the  propellant  charge  in  the  main  igniter  chamber,  it  is  diffi- 
cult to  control  the  surface  area  of  the  shavings  and  thus  hard  to  con- 
trol the  igniter  pressure  pulse.  In  order  to  achieve  high  reproduci- 
bility of  the  pressurization  process  in  the  crack  combustion  chamber, 
the  propellant  shavings  were  deleted  from  the  booster  chamber.  To 
facilitate  and  ensure  a successful  ignition  of  the  propellant  charge  in 
the  main  igniter  chamber,  the  propellant  charge  was  coated  with  an 
igniter  paste  (T1-B-KC104]  suggested  by  Mr.  H.  B.  Mathes  of  NWC.*  This 
set-up  has  proven  to  be  successful  to  generate  reproducible  pressuriza- 
tion rates,  on  the  order  of  10^  atm/s  at  the  opening  of  the  crack.  The 
initial  pressurization  rate  in  the  crack  combustion  chamber  can  be 
controlled  by  varying  the  igniter  charge  weight,  the  area  of  the  inlet 
perforated  nozzle,  and  the  area  of  the  exit  nozzle  in  the  combustion 
chamber.  The  exploded  view  of  the  present  igniter  system  is  shown  in 
Figure  2.  The  complete  assembly  view  of  the  igniter  system  mounted 
together  with  the  crack  combustion  chamber  is  shown  in  Figure  3. 


A - Solenoid  armature 
B - Percussion  element 
C - Retainer  for  percus- 
sion element 


D - Booster  chamber 
E - Propellant  igniter  charge 
F - Igniter  chamber 
G - Perforated  nozzle  plate 


FIGURE  2.  Exploded  View  of  Solid  Propellant  Igniter  System. 


* 

Personal  conmunicatlon  with  H.  B.  Mathes,  NWC,  19  July  1977. 
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FIGURE  3.  Assembly  View  of  the  Igniter  and  Crack  Combustion  Chamber. 


SOLID  PROPELLANT  CRACK  SAMPLES 

The  propellant  samples  used  in  earlier  firings  were  cast  from 
smokeless  propellant  supplied  by  NWC.  The  method  of  casting  and  the 
arrangement  of  casting  mold  with  teflon-coated  mandrel  were  discussed  in 
the  previous  annual  report.1  Since  the  smokeless  propellant  is  very 
viscous  in  its  liquid  state,  it  is  difficult  to  eliminate  bubbles  from 
the  cast  sample.  Therefore,  a rubber-based  composite  propellant,  des- 
ignated as  Propellant  A,  has  been  used  to  replace  the  smokeless  propel- 
lant in  our  experimental  investigations.  Two  different  mandrel  sizes 
have  been  used  in  the  casting  process;  one  to  form  a large  crack  with  an 
average  gap  width  of  0.52  cm  and  a length  of  20  cm,  and  the  other  to 
form  a small  crack  with  an  average  gap  width  of  0. 10  cm  and  a length  of 
15  cm.  The  case  sample  for  a small  crack  is  shown  in  Figure  4. 
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FIGURE  4.  Cast  Propellant  Crack  Sample. 


Test  firings  were  conducted  with  the  cast  samples  of  both  large  and 
small  cracks.  In  these  firings,  propellant  debond  at  the  interface  of 
the  sample  and  the  brass  mold  was  observed.  The  flame  propagation  along 
the  debond  is  given  in  detail  in  Discussion  of  Results.  The  cause  for 
the  debond  is  believed  to  be  due  to  the  bonding  weakness  at  the  inter- 
face of  the  two  different  materials.  To  reduce  the  possibility  of 
debond,  the  propellant  specimen  was  trimmed  off  at  the  crack  opening  and 
the  interface  corners  of  the  recessed  cavity  were  reinforced  with  epoxy. 
This  procedure  has  proven  to  be  successful  in  eliminating  the  debond 
problem.  As  a result  of  recessing  the  entrance  plane  of  the  crack 
sample,  the  length  of  the  crack  is  slightly  less  than  the  initially 
designed  mandrel  length. 


There  are  several  undesirable  features  of  the  cast  propellant 
samples: 

1.  The  gap  width  of  the  crack  is  usually  nonuniform  along  the 
crack  length,  even  though  the  mandrel  is  of  uniform  thickness.  The 
percentage  of  gap  width  variation  could  be  as  high  as  60%  for  small 
cracks.  It  is  very  unlikely  to  have  crack  samples  of  same  initial 
geometrical  configurations;  this  presents  problems  in  repeating  a given 
experiment. 

2.  The  mandrels  may  deflect  from  the  centerline  during  the  cast- 
ing process.  This  deflection  results  in  the  crack  gap  not  aligning  with 
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the  pressure  gauges.  Some  cast  samples  have  crack  gaps  that  completely 
miss  the  pressure  gauge  port  near  the  crack  tip  location. 

3.  A cast  propellant  sample  may  have  recessed  areas  on  its  sur- 
face; if  these  areas  are  not  filled  with  epoxy,  a flame  or  hot  gases  may 
penetrate  into  these  defect  areas. 

4.  The  cast  propellant  crack  has  fuel-rich  surfaces  which  differ 
from  the  crack  surfaces  generated  by  the  mechanical  damage  of  the  pro- 
pel lant . 

In  view  of  these  undesirable  features  of  the  cast  cracks,  a pro- 
cedure has  been  developed  to  fabricate  propellant  crack  samples  with  a 
slit-type  crack.  In  this  method,  the  brass  mold  is  replaced  by  the 
assembly  of  two  brass  legs  and  an  upper  metal  insert.  The  crack  is 
formed  in  between  two  slabs  of  propellant.  The  procedure  to  prepare  a 
sample  is  performed  by  first  adhesively  bonding  an  oversized  propellant 
slab  to  each  leg  of  the  modified  brass  mold.  Next,  each  leg  of  the  mold 
is  mounted  on  top  of  a machinist's  indexing  table  and  secured  by  a 
retaining  system.  The  indexing  table  is  attached  to  a drill  press  and 
all  exposed  surfaces  of  the  propellant  slab  are  precisely  cut  to  t. 
desired  dimensions  by  an  end  mill  (see  Figure  5) . A low-pressure  au 
line  is  used  during  the  operation  to  clear  propellant  shavings.  The 
shavings  are  collected  and  later  burnt.  The  operator  of  the  drill  press 
wears  a proximity  fire  suit  as  a safety  precaution.  The  final  crack 
configuration  is  shown  in  Figure  6. 

The  slit-type  propellant  cracks  have  the  advantage  of  better  con- 
trolled gap  width  distributions,  thus  are  more  suitable  for  obtaining  a 
high  degree  of  repeatability.  Another  advantage  of  this  crack  arrange- 
ment is  the  ability  to  form  extremely  small  gap  widths.  Also,  the 
length  of  the  crack  can  be  easily  varied  by  changing  the  length  of  the 
upper  insert  piece.  The  machined  propellant  surface  is  closer  to  the 
crack  surface  generated  by  mechanical  damage  of  the  propellant  and, 
therefore,  the  surface  gasification  temperature  is  that  of  the  propel- 
lant rather  than  the  fuel-rich  surface. 

The  slit-type  propellant  cracks  also  have  the  flexibility  to  study 
the  crack  tip  ignition  phenomena.  This  can  be  achieved  by  adhesively 
bonding  a piece  of  propellant  to  the  bottom  surface  of  the  upper  insert 
piece.  A propellant  crack  sample  of  this  modified  configuration  is 
shown  in  Figure  7. 
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FIGURE  5.  Milling  of  Solid  Propellant  Slabs  to  Form 
Slit-Type  Propellant  Cracks. 
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FIGURE  6.  Slit-Type  Propellant  Crack  Sample 


FIGURE  7.  Slit-Type  Propellant  Crack  Sample  Modified  to 
Study  Crack  Tip  Ignition  Phenomena. 
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In  order  to  check  the  one  dimensionality  of  the  flame  spreading 
process  in  narrow  cracks  of  various  crack  width-to-thickness  ratios,  a 
new  arrangement  of  the  crack  configuration  has  been  designed.  This  new 
configuration  is  shown  in  Figure  8.  The  crack  is  formed  between  a 
propellant  slab  slightly  regressed  below  the  leg  assembly  and  the  inner 
surface  of  the  plexiglas  sacrificial  window  which  is  mounted  on  top  of 
the  leg  assembly.  This  type  of  arrangement  allows  the  observation  of 
the  flame  spreading  process  from  the  front  view  of  the  crack  surface  and 
thus  determines  the  uniformity  of  the  propagating  flame  front.  Two 
pressure-time  traces  will  be  recorded:  one  at  the  opening  of  the  crack 
and  one  near  the  tip  of  the  crack  through  a slot  in  the  upper  insert 
piece.  Although  the  fabrication  of  the  components  for  this  mold  assem- 
bly has  been  completed,  the  experimental  firings  were  not  performed 
since  they  were  not  included  in  the  scope  of  this  reporting  period. 
However,  tests  will  be  conducted  in  the  immediate  future  under  the 
continuing  contract. 


FIGURE  8.  Propellant  Crack  Sample  for  the 
Study  of  One-Dimensionality  of  the  Flame 
Propagation  Phenomena. 

m 
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CAMERA  SYSTEM 

A new  camera  has  been  incorporated  into  the  data  acquisition  system 
to  provide  for  time-correlation  between  the  pressure  traces  and  flame- 
spreading data.  The  camera  (100'  HYCAM)  and  associated  LED  drive  unit 
are  on  loan  from  NSWC/DL.  The  event  synchronizer  circuit  of  the  camera 
was  used  to  trigger  the  solenoid-actuated  firing  pin  in  the  igniter. 

This  special  circuit  was  found  to  be  very  advantageous  for  capturing  the 
flame  spread  transient  and  ensuring  that  the  camera  was  operating  at  the 
preselected  framing  rate  prior  to  ignition.  Unfortunately,  the  100' 
HYCAM  does  not  operate  at  fast  enough  frame  rates  (3000  pps  maximum  with 
regulated  speed  control)  to  give  sufficient  time  resolution  of  the  high- 
speed flame  spread  interval.  This  deficiency  has  been  recently  allevi- 
ated with  a 400'  HYCAM  camera  acquired  through  an  interdepartmental  loan 
at  Pennsylvania  State  University.  This  camera,  which  has  features 
similar  to  those  of  the  100'  HYCAM,  can  achieve  frame  rates  up  to  10,000 
pps  (5,000  pps  for  speed  regulated  mode).  The  full-frame  and  accessory 
half-frame  heads  of  the  camera  have  been  modified  to  interface  with  a 
high-precision,  externally-driven,  LED  timing  light  generator.  To 
facilitate  future  test  firings,  a new  LED  driver  unit  similar  to  the  one 
on  loan  from  NSWC/DL  has  been  designed  and  fabricated.  Figure  9 shows 
an  overall  view  of  the  filming  process. 


TEST  PROCEDURE 

After  the  propellant  crack  sample  has  been  prepared  and  mounted 
into  the  crack  combustion  chamber,  the  crack  geometry  is  carefully 
measured  and  recorded  on  a data  sheet.  The  crack  combustion  chamber  is 
then  assembled  and  installed  in  the  test  cell.  The  piezoelectric  pres- 
sure transducers  are  then  mounted  into  the  crack  combustion  chamber  and 
cooling  water  lines  are  connected  to  the  transducer  water-cooled 
adapters. 

To  provide  maximum  safety  precautions  and  to  minimize  the  possi- 
blity  of  errors,  the  remainder  of  the  firing  is  performed  by  following 
the  checklist  given  in  Appendix  A.  As  a result  of  this  procedure,  a 
high  degree  of  successful  firings  have  been  obtained.  A sample  data 
sheet  used  in  the  crack  combustion  experiments  is  shown  in  Appendix  B. 


DISCUSSION  OF  RESULTS 


PRESENTATION  AND  DISCUSSION  OF  EXPERIMENTAL  DATA 

During  this  reporting  period,  a number  of  test  firings  have  been 
conducted  to  determine  the  influence  of  both  the  solid  propellant  ig- 
niter and  the  crack  geometry  on  the  flame-spreading  and  pressurization 
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in  propellant  crack  samples.  The  propellant  samples  were  either  cast  by 
NWC  or  machined  to  a slit-type  crack  at  Pennsylvania  State  University. 
The  test  runs  discussed  in  this  section  were  all  conducted  by  using  the 
recently  developed  solid  propellant  igniter  system.  The  strength  of  the 
igniter  was  varied  by  changing  the  propellant  charge  weight  in  the 
igniter  chamber. 


In  the  earlier  tests,  propellant  debond  at  its  interface  with  the 
brass  mold  was  observed.  DNC  Test  No.  12  represents  a typical  case  of 
the  unplanned  debond  experiments.  In  this  test  the  crack  sample  was 
cast  from  smokeless  propellant.  The  gap  width  was  uniform  at  a value  of 
0.51  cm  and  the  length  of  the  crack  was  20  cm.  The  igniter  system  was 
composed  of  an  FA-34  primer,  AF  composite  propellant  shavings  as  booster, 
and  a rectangular  charge  of  AF  propellant  with  a weight  of  2.284  g. 


FIGURE  9.  Overall  View  of  the  Filming  Process. 
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Figure  10  shows  the  set  of  pressure-time  traces  recorded  at  various 
locations  along  the  crack.  Time-correlation  of  the  pressure  and  flame 
spreading  data  is  achieved  with  the  common-time  pulse.  The' pulse  was 
generated  by  an  LED  driver  when  the  igniter  solenoid  was  energized.  As 
measured  by  Gauge  2 (G2),  the  igniter  produced  an  initial  peak  pressure 
of  120  atm  at  a rate  of  1.76  x 10^  atm/s  in  the  main  chamber  of  the  test 
rig.  The  shape  of  the  p-t  trace  near  the  first  peak  of  G2  closely 
matches  the  profiles  obtained  from  igniter  tests  with  a simulated  inert 
propellant  sample.  This  implies  that  the  initial  pressure  pulse  is 
solely  due  to  the  igniter  output.  During  this  initial  pressure  pulse 
interval,  the  pressurization  rates  differ  from  gauge-to-gauge.  All 
gauge  pressures  rise,  due  to  gasification  processes,  as  the  combustion 
process  continues  until  a second  peak  is  reached. 

The  high-speed  film  (2540  pps)  revealed  that  ignition  occurred  at 
the  opening  of  the  debond  at  approximately  1.2  ms  after  the  first  glimpse 
of  hot  igniter  gas  in  the  chamber.  The  ignition  front  then  proceeded 
rapidly  along  a developing  debond  between  the  propellant  and  the  brass 
mold.  Flame  spreading  in  the  crack  occurs  at  a later  time.  Certain 
frames  from  this  film  were  selected  to  show  important  events  observed  in 
the  combustion  process  (see  Figure  11) . Observations  and  interpreta- 
tions of  the  combustion  phenomena  associated  with  each  picture  are 
listed  in  Table  1. 
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FIGURE  10.  Experimental  Pressure-Time  Traces 
from  DNC  Test  No.  12. 
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FIGURE  11.  Motion  Pictures  of  the  Ignition  and  Combustion 
Phenomena  in  a Solid  Propellant  Crack  with  Unplanned  Debond 
in  DNC  Test  No.  12. 
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FIGURE  11.  (Contd.) 
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TABLE  1.  Physical  Interpretation  of  the  Film 
of  DNC  Test  No.  12. 


Picture  No. 

Time,  ms 

Observations  and  interpretations 

1 

0.0 

Hot  igniter  gas  enters  the  rectangular  duct  in 
the  bottom  of  the  main  chamber  and  begins  to 
penetrate  into  the  crack 

2 

0.39 

Brightness  in  the  debond  region  showing  pene- 
tration of  hot  gases  into  the  debond  crack 

3 

0.79 

The  developing  debond  between  the  propellant 
and  the  brass  mold  closes  the  crack  opening 

4 

1.18 

Possible  ignition  near  the  opening  of  the 
debond 

5 

1.57 

There  is  a propagation  of  debond  crack 

7 

2.36 

Gas  pressure  starts  to  open  the  crack  and 
close  the  debond 

9 

3.15 

Extinction  occurs  in  a part  of  the  debond  due 
to  closure  of  that  portion  of  the  debond 

11 

3.94 

Severe  burning  in  the  debond  once  again  closes 
the  crack  opening 

13 

4.72 

Severe  burning  continues  in  the  debond 

15 

5.51 

The  crack  entrance  begins  to  open 

18 

6.69 

The  crack  is  closed 

21 

7.87 

The  crack  is  open,  yet  there  is  still  burning  ii 
the  debond 

22 

8.27 

The  crack  is  closed 

23 

8.66 

The  crack  is  open 

25 

9.45 

The  crack  is  open,  yet  there  is  still  burning  ii 
the  debond 

27 

10.24 

There  is  severe  burning  in  the  crack.  Gas  start 
to  flow  into  the  gap  between  the  propellant 
sample  and  the  chamber  window 

49 

18.90 

The  debond  piece  becomes  thin  due  to  burning  on 
both  surfaces 

31 


NWC  TP  6049 


The  next  test  firing  data  selected  for  discussion  was  DNC  Test  No. 

16.  The  crack  sample  used  was  cast  by  NWC  from  Type  A propellant.  The 
average  gap  width  was  0.14  cm  and  the  length  of  the  crack  was  14.24  cm. 

The  igniter  was  the  same  as  DNC  Test  12  except  the  propellant  charge 

weight  for  this  test  firing  was  1.925  g.  The  debond  problem  was  elim-  1 

inated  by  utilizing  the  method  described  in  the  section  entitled 
"Experimental  Set-Up  of  the  Test  Rig." 

Figure  12  shows  the  measured  pressure-time  traces  along  the  crack. 

The  average  pressurization  rate  at  G2  before  reaching  the  peak  pressure 
is  approximately  0.7  x 10^  atm/s.  It  is  interesting  to  note  that  at 
this  high  pressurization  rate,  the  pressure  irregularities  at  downstream 
gauge  locations  immediately  follow  the  first  discernible  pressure  rate 
at  G2.  These  irregularities  are  believed  to  be  caused  by  strong  wave 
propagation  and  reflection  phenomena  inside  the  crack.  These  pressure 

irregularities  continue  to  appear  until  the  pressure  peak  region  is  < 

over.  It  is  also  important  to  note  that  the  pressure  trace  at  each 
gauge  location  displays  only  one  pressure  peak  region. 

The  pictures  showing  the  important  gas  penetration  and  flame  spread- 
ing event  are  given  in  Figure  13.  The  detailed  observations  of  these 
pictures  are  listed  in  Table  2.  Generally  speaking,  the  flame  spreading 
information  obtained  in  this  firing  seems  to  suggest  that  when  the  com- 
pression wave  generated  by  the  igniter  system  is  too  strong,  the  region 
near  the  crack  tip  may  be  ignited  earlier  than  the  region  closer  to  the 
entrance  of  the  crack.  This  phenomena  produces  flame  spreading  in  the 
reversed  direction  near  the  crack  tip  which  in  turn  can  lead  to  the 
convergence  of  two  propagating  flame  fronts. 
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FIGURE  12.  Experimental  Pressure-Time  Traces 
from  DNC  Test  No.  16. 
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FIGURE  13.  (Contd.) 


TABLE  2.  Physical  Interpretation  of  the  Film 
of  DNC  Test  No.  16. 


Picture  No. 

Time,  ms 

Observations  and  Interpretations 

1 

0.0 

First  glimpse  of  igniter  gas  at  the  crack 
entrance 

2 

0.38 

Hot  gases  start  to  penetrate  the  crack  near  the 
entrance 

3 

0.77 

Hot  gases  are  noticed  in  the  downstream  portion 
of  the  crack 

4 

1.15 

Possible  ignition  in  the  region  between  9 to 

13  cm  from  the  crack  opening 

5 

1.54 

Flame  propagation  in  the  reversed  direction 
toward  the  crack  opening 

6 

1.92 

Severe  burning  near  the  center  of  the  crack 

7 

2.30 

The  area  at  the  crack  entrance  is  substantially 
reduced 

8 

2.69 

Severe  burning  inside  the  crack 

15 

5.38 

After  severe  burning  has  continued  for  2 ms, 
the  crack  entrance  is  enlarged  again 
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The  compression  wave  generated  ignition  involves  a new  mechanism  of 
ignition  other  than  the  conventional  ignition  caused  by  convective 
heating  of  the  propellant  surface.  Currently,  this  new  mechanism  has 
not  been  considered  in  the  crack  combustion  model;  therefore,  the  con- 
vergent flame  propagation  phenomena  cannot  be  predicted  by  our  computer 
code.  This  new  mechanism  of  ignition  in  the  authors’  opinion  deserves 
future  investigations. 

In  DNC  Test  No.  18  the  propellant  sample  was  again  cast  by  NWC  from 
Type  A propellant.  The  average  gap  width  was  0.16  cm  and  the  crack 
length  was  14.73  cm.  The  same  igniter  was  used  but  with  a propellant 
charge  weight  of  1.966  g.  The  pressure-time  traces  are  shown  in  Fig- 
ure 14.  The  pressurization  rate  of  G2  at  the  onset  of  pressurization  is 
approximately  2.2  x 10^  atm/s  which  is  nearly  three  times  the  rate  of 
DNC  Test  No.  16.  Again,  trhe  strong  pressure  irregularities  following 
the  first  discernible  pressure  rise  are  noticed  in  the  downstream 
gauges.  There  are  mild  second  peak  pressures  reached  at  time  consid- 
erably longer  than  initial  pressurization  time. 

The  pictures  describing  the  observed  gas  penetration,  ignition,  and 
flame  propagation  phenomena  are  given  in  Figure  15.  The  interpretations 
' of  these  pictures  are  given  in  Table  3.  Similar  to  DNC  Test  No.  16,  the 
tip  region  seems  to  be  ignited  before  the  flame  propagates  from  the 
crack  entrance  to  the  tip,  resulting  in  flame  propagation  in  both  direc- 
tions . 
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FIGURE  14.  Experimental  Pressure-Time  Traces 
from  DNC  Test  No.  18. 
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FIGURE  15.  Motion  Pictures  of  the  Ignition  and 
Combustion  Phenomena  in  a Solid  Propellant  Crack 
of  DNC  Test  No.  18. 
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TABLE  3.  Physical  Interpretation  of  the  Film 
of  DNC  Test  No.  18 


Picture  No. 

Time,  ms 

Observations  and  interpretations 

1 

0.0 

Hot  igniter  gases  enter  the  combustion  chamber 
and  start  to  penetrate  the  crack 

2 

0.38 

Hot  gases  are  noticed  near  the  crack  tip  indi- 
cating possible  ignition  in  this  region. 

3 

0.77 

Flame  propagates  in  both  directions 

6 

1.92 

Severe  burning  near  the  region  3.5  cm  from  the 
crack  opening 

9 

3.08 

The  crack  tip  region  and  the  3.5  cm  region  are 
burning  most  severely,  while  other  regions 
remain  less  intense 

DNC  Test  No.  22  was  conducted  by  using  the  slit-type  crack  sample 
fabricated  from  Propellant  A.  The  gap  width  of  the  crack  was  a uniform 
value  of  0.19  cm  and  the  length  was  14  cm.  The  igniter  consisted  of  an 
FA-34  primer  and  an  igniter  charge  of  Propellant  A at  a weight  of  2.066 
g which  was  coated  with  Ti-B-KClO^  ignition  paste.  Since  the  igniter 
charge  was  made  of  Propellant  A rather  than  the  more  energetic  AF  com- 
posite propellant  and  also  there  were  no  propellant  shavings  in  the 
booster  chamber,  the  initial  chamber  pressurization  rate  was  much  lower 
than  the  cases  discussed  earlier.  The  initial  pressurization  rate  for 
this  test  was  0.166  x 10-*  atm/s. 

The  pressure-time  traces  at  the  crack  opening  and  the  crack  tip  for 
this  test  firing  are  shown  in  Figure  16.  Due  to  the  fact  that  the  crack 
gap  width  is  larger  than  those  in  DNC  Test  No.  16  and  18,  and  also  the 
initial  pressurization  rate  for  this  test  is  much  lower,  the  downstream 
gauge  does  not  have  the  pressure  irregularities  recorded  in  the  two 
previously  discussed  tests.  The  pressure  gradient  along  the  crack  is 
very  small  even  in  the  initial  pressurization  interval;  this  is  con- 
sistent with  the  theoretical  findings  for  cracks  of  relatively  large  gap 
widths  ignited  with  relatively  low  chamber  pressurization  rates. 

The  pictures  showing  the  important  events  occurring  in  the  time 
interval  of  interest  are  shown  in  Figure  17.  The  interpretations  of  the 
observed  phenomena  for  these  pictures  are  listed  in  Table  4.  The  ad- 
vancement of  the  flame  front  can  be  seen  clearly  from  these  pictures  and 
the  flame  deceleration  near  the  closed  end  of  the  crack  agrees  with  the 
observations  made  separately  by  Belyaev,  et  al . and  Taylor.14*12 

12J.  W.  Taylor.  "The  Burning  of  Secondary  Explosive  Powders  by  a 
Convective  Mechanism,"  Trans.  Farad.  Soc.  , Vol.  58  (1962),  p.  561. 
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FIGURE  16.  Experimental  Pressure-Time 
Traces  from  DNC  Test  No.  22. 


COMPARISON  OF  THEORY  WITH  EXPERIMENT 

Figure  18  shows  a typical  set  of  experimental  and  theoretical 
p-t  traces.  The  experimental  p-t  traces  shown  were  measured  from  DNC 
Test  No.  21.  The  crack  specimen  for  this  test  firing  was  made  of  pro- 
pellant Type  A cast  by  NWC.  The  length  of  the  crack  was  14.6  cm  and  the 
gap  width  at  the  crack  opening  was  0.089  cm,  increasing  gradually  to 
0.140  cm  at  the  tip.  The  theoretical  p-t  traces  shown  were  obtained 
from  the  crack  combustion  code  using  the  input  data  listed  in  Table  5. 
Time  zero  in  Figure  18  is  defined  as  the  time  corresponding  to  the  first 
glimpse  of  igniter  gas  on  the  high-speed  film. 

Several  special  features  were  observed  by  comparing  the  experi- 
mental p-t  traces  measured  at  various  locations.  The  first  discernible 
pressure  rise  in  G2  occurs  very  shortly  (At  <0.2  ms)  after  the  first 
glimpse  of  hot  igniter  gas.  Each  gauge  measures  two  peaks;  the  first 
peak  is  essentially  caused  by  the  discharge  of  the  igniter  material  and 
the  second  peak  is  caused  by  the  burning  inside  the  propellant  crack. 
Although  the  pressures  at  the  first  peak  do  not  differ  significantly  for 
the  three  gauge  locations,  the  rates  of  pressurization  to  achieve  the 
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FIGURE  17.  Motion  Pictures  of  the  Ignition 
and  Combustion  Phenomena  in  a Solid  Propel- 
lant Crack  of  DNC  Test  No.  22. 
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FIGURE  17.  (Contd.) 


TABLE  4.  Physical  Interpretation  of  the  Film 
of  DNC  Test  No.  22. 


Picture  No.  I Time,  ms 


Observations  and  interpretations 


0.0  Hot  igniter  gas  enters  the  rectangular  duct  in 
the  bottom  of  the  main  chamber 

0.2,0. 4 Hot  igniter  gas  penetrates  into  the  crack 


0.6,0. 8 


1.0, 1.2 


1.8, 2.0 


2. 6, 2. 8 


Brightness  is  observed  near  the  crack  tip.  The 
entrance  portion  of  the  crack  may  be  ignited 

Hot  igniter  gases  penetrate  toward  the  crack  tip. 
The  flame  front  advances  from  5 cm  to  6.5  cm 
measured  from  the  crack  opening 

The  flame  front  advances  from  7.5  to  8.5  cm 
measured  from  the  crack  opening 

Flame  propagation  continues.  The  flame  front 
advances  from  10  to  11  cm  measured  from  the 
crack  opening 

Brightness  near  the  crack  tip  suggests  ignition 
in  this  region.  Flame  front  is  decelerated  near 
the  crack  tip  region  but  still  advances  from 
11.5  to  12  cm 

High  intensity  gas  fills  the  entire  length  of 
the  crack 
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TABLE  5.  Computer  Program  Input  Variables 
in  the  Simulation  of  DNC  Test  No.  21. 
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-S 

0.0  cm 

p 

XL 

x 

14.61  cm 

B 

X 

- 

0.0  gf/g 

b 

X 

1.0  cm'Vg 

M 

w 

X 

26.10  g^/g-mole 

Y 

“ 

1.21 

P„_ 

3 

1.71  g/cm^ 

pr 

A 

3 

0.80  x 10  ^ cal/s-cm-°K 

Pr 

2 

a 

X 

0.0018  cm  /s 

pr 

es/dh 

X 

0.025 

Tf 

X 

3000  K 

T , 

X 

860  K 

cri 

T . 

X 

920  K 

ign 

T . 

X 

298.5  K 

i 

a 

X 

0.9359  x 10'2(— )/(^x)n 

s i 

n 

X 

0.4108  Cm 

K 

e 

X 

0.70  cm^-K/cal 

6 

X 

53 

T 

X 

2400  K 

c 

first  peak  pressure  do  differ  significantly.  G4  has  the  highest  rate  of 
pressurization.  The  first  discernible  pressure  rise  occurs  consecutive- 
ly at  the  G2,  G3,  and  G4  locations.  After  the  first  peak,  pressure 
decays  at  all  three  gauge  locations.  There  are  more  pressure  irregu- 
larities at  G4  than  the  other  two  gauges  in  the  time  period  immediately 
following  the  first  peak.  The  pressure  decay  after  the  first  peak  is 
believed  to  be  caused  by  a net  decrease  of  gaseous  mass  inside  the 
crack.  This  mass  reduction  results  from  a higher  rate  of  gas  flowing 
out  of  the  crack  than  the  rate  of  gasification  of  propellant  inside  the 
crack.  The  pressure  differences  between  various  locations  along  the 
crack  become  smaller  and  smaller  as  the  second  peak  is  approached. 

After  the  second  peak,  the  pressure  decays  continuously  due  to  the 
burnout  of  the  propellant  sample. 
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The  theoretical  p-t  traces,  shown  as  dashed  lines  in  Figure  18,  are 
in  good  qualitative  agreement  with  the  experimental  results.  The  time 
corresponding  to  the  first  discernible  pressure  rise  at  G2,  G3,  and  G4 
increases  consecutively.  Similar  to  the  experimental  traces,  the  cal- 
culated initial  pressurization  rates  differ  significantly  from  each 
other  at  the  three  gauge  locations.  The  pressurization  rate  at  G4  is 
much  higher  than  the  rate  of  pressurization  at  the  opening  of  the  crack. 
Also,  the  calculated  p-t  traces  display  pressure  irregularities  shortly 
after  the  first  peak  near  the  crack  tip.  These  spikes  are  due  to  compli- 
cated compression  wave  phenomena  near  the  crack  tip  region.  Furthermore, 
in  the  same  manner  as  the  experimental  pressure  traces,  the  theoretically 
calculated  pressure  distribution  along  the  crack  becomes  more  uniform  as 
the  second  peak  is  approached. 

In  spite  of  the  qualitative  similarities  discussed  above,  there  is 
some  quantitative  disagreement  between  theory  and  experiment.  The 
calculated  pressures  at  G3  and  G4  are  higher  than  the  experimentally 
measured  pressure  for  a large  portion  of  the  time  period  being  examined. 
The  theoretically  calculated  first  discernible  pressure  rise  at  G4 
occurs  before  the  experimentally  observed  time.  Also,  the  calculated 
first  peaks  for  G3  and  G4  are  greater  than  the  experimentally  measured 
values. 


The  quantitative  disagreements  between  theory  and  experiment  can  be 
caused  by  several  factors  which  are  not  considered  in  the  present  model. 
The  current  version  of  the  crack  combustion  code  does  not  account  for 
any  mechanical  deformation  of  the  crack.  The  initial  pressure  pulse 
from  the  igniter  into  the  rectangular  duct  in  the  bottom  of  the  main 
combustion  chamber  could  mechanically  push  the  bottom  surface  of  the 
propellant  specimen  inward  and  thereby  decrease  the  crack  gap  width  near 
the  crack  opening.  Based  upon  an  incompressible  assumption,  mechanical 
deformation  of  the  solid  propellant  can  reduce  the  crack  gap  width  by 
50%.  This  decrease  in  flow  area  would  result  in  increased  flow  resis- 
tance into  the  crack.  The  flow  rate  of  penetrating  gases  is  reduced  by 
both  a decrease  in  flow  area  and  a decrease  in  discharge  coefficient. 

The  change  of  propellant  contour  near  the  crack  opening  may  substan- 
tially reduce  the  discharge  coefficient.  If  the  crack  combustion  code 
considered  the  mechanical  deformation  of  the  propellant,  then  the  first 
discernible  pressure  rise  at  the  tip  would  occur  at  a later  time.  Also, 
if  the  compressibility  of  the  solid  propellant  was  considered  in  the 
model,  a lower  pressure  would  be  obtained  in  the  crack.  This  indicates 
the  necessity  of  coupling  the  crack  combustion  code  with  a structure 
analysis  program. 

In  general,  the  theoretical  results  obtained  from  the  crack  combus- 
tion code  are  in  good  qualitative  agreement  with  experimental  observa- 
tion. Although  the  detailed  quantitative  comparison  of  the  theoretical 
results  with  experimental  data  showed  noticeable  differences,  the  cal- 
culated pressure-time  traces  and  the  pressurization  rates  after  the 
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first  discernible  pressure  rise  are  in  reasonable  agreement  with  the 
data.  It  is  our  belief  that  quantitative  match  would  be  improved  after 
coupling  the  crack  combustion  code  with  the  structure  analysis  program 
(NONSAP) . 

Figure  19  describes  .the  pressure  distribution  along  the  crack  at 
various  times.  Due  to  the  initial  pressurization  caused  by  the  igniter, 
the  pressure  at  the  crack  opening  is  greater  than  that  near  the  crack 
tip.  The  pressure  gradients  at  times  less  than  1 ms  are  favorable  for 
the  penetration  of  hot  gases  into  the  propellant  crack.  Between  1 ms 
and  1.25  ms,  there  is  a substantial  increase  in  pressure  in  a large 
portion  of  the  crack.  This  increase  is  caused  by  the  gasification  of 
propellant  inside  the  crack.  The  accumulation  of  gases  inside  the  crack 
produces  a reverse  pressure  gradient  near  the  crack  opening.  The  pro- 
duct gases  are  driven  out  of  the  crack  by  the  reverse  pressure  gradient. 
The  outflow  of  product  gases  causes  a decrease  in  pressure  inside  the 
crack.  This  is  why  there  are  valleys  on  the  p-t  traces  shown  in  Figure 
18.  At  the  same  time,  the  outflow  of  product  gases  greatly  decelerates 
the  flame  spreading  process  at  the  center  of  the  crack. 

The  calculated  locus  of  ignition  front  is  superimposed  as  a dashed 
line  in  Figure  19.  The  locus  of  ignition  front  indicates  rapid  flame 
spreading  from  the  crack  opening  to  7.7  cm.  The  flame  spreading  abruptly 
decelerates  near  the  center  of  the  crack  due  to  the  reason  given  above. 
Additional  explanation  of  the  flame  deceleration  is  given  below  in  the 
discussion  of  the  calculated  velocity  distributions. 

The  flame  spreading  and  combustion  phenomena  for  DNC  Test  No.  21 
were  recorded  with  a 16-mm  high-speed  motion  picture  camera  at  an  average 
rate  of  5100  pps.  A number  of  frames  from  this  film  were  selected  to 
show  the  important  events  in  the  overall  combustion  process  (see  Fig- 
ure 20).  A description  of  the  observed  phenomena  is  given  in  Table  6. 
Similar  to  the  theoretical  results,  the  ignition  front  advances  at  a 
high  speed  through  the  first'  half  of  the  crack  length  and  then  the  flame 
front  becomes  stationary  near  the  center  of  the  crack. 

Figure  21  shows  the  calculated  gas  velocity  profiles  that  exist  in 
the  crack  at  various  times.  At  the  time  of  0.5  ms  the  hot  gases  are 
flowing  into  the  crack  at  a high  velocity  because  of  a favorable  initial 
pressure  gradient.  At  1 ms  the  velocity  of  the  inflowing  gases  de- 
creases as  the  pressure  gradient  at  the  crack  opening  decreases.  As  the 
flame  front  progresses  along  the  crack,  the  resultant  gasification 
processes  cause  the  pressure  gradient  at  the  entrance  of  the  crack  to 
reverse.  At  2 ms,  the  product  gases  are  driven  out  of  the  crack  by  the 
reversed  pressure  gradient.  The  decrease  in  velocity  of  the  outflowing 
gases  between  6 and  10.5  ms  is  probably  due  to  the  increase  in  port  area 
at  the  crack  entrance  which  will  be  discussed  later.  During  the  entire 
pressurization  process,  the  gas  velocity  near  the  closed  end  of  the 
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2.73  ms  5.07  ms  6.24  ms  8.97  ms  12.87  ms  15.21  ms 


FIGURE  20.  Motion  Pictures  of  the  Ignition  and  Combustion 
Phenomena  In  a Solid  Propellant  Crack  of  DNC  Test  No.  21. 
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TABLE  6.  Physical  Interpretation  of  the  Film 
of  DNC  Test  No.  21. 


Picture  No. 

Time,  ms 

Observations  and  interpretations 

1 

0.0 

Hot  igniter  gas  enters  the  rectangular  duct 
in  the  bottom  of  the  main  chamber 

3 

0.39 

The  hot  igniter  gas  penetrates  into  the 
crack 

4 

0.59 

High  intensity  gases  appear  at  the  tip  of 
the  crack 

5 

0.78 

High  intensity  gases  disappeared  from  the 
tip 

9 

1.56 

The  booster  piece,  attached  to  the  left  leg 
at  the  crack  entrance,  burns  severely 

15 

2.73 

High  intensity  light  in  the  region  near  the 
crack  opening 

27 

5.07 

After  a period  of  decrease  in  light  intensil 
the  intensity  increases  again 

32 

6.24 

The  flame  is  not  propagating 

47 

8.97 

The  flame  still  has  not  advanced  further  up 
the  crack 

67 

12.87 

The  flame  is  not  propagating  up  the  crack 

79 

15.21 

The  entire  crack  is  ignited 

crack  is  very  low;  this  is  mainly  due  to  the  stagnation  point  condition 
at  the  crack  tip.  Another  reason  for  the  low  velocity  in  this  region  is 
the  lack  of  a favorable  pressure  gradient  as  the  driving  potential.  The 
low  flow  velocity  significantly  reduces  the  convective  heating  of  the 
propellant  surface.  Hence,  a much  longer  time  is  required  for  the 
propellant  surface  near  the  crack  tip  to  achieve  a state  of  ignition. 
This  delay  is  readily  evident  in  both  the  calculated  locus  of  ignition 
shown  in  Figure  19  and  the  high-speed  film  in  Figure  20. 
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FICURE  21.  Calculated  Velocity  Distribution  at 
Various  Times  (Simulation  of  DNC  Test  No.  21). 
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Figure  22  gives  the  calculated  gas  temperature  distribution  along 
the  crack  at  various  times.  The  temperature  gradients  at  0.1  and  0.5  ms 
are  very  pronounced  near  the  opening  of  the  crack  indicating  the  pene- 
tration of  hot  gases  into  the  propellant  crack.  Figure  22  shows  that 
the  gas  temperature  is  lowest  near  the  tip  of  the  crack.  The  local  gas 
temperature  is  low  because  no  favorable  pressure  gradient  exists  to 
cause  the  hot  gases  to  flow  into  the  stagnation  region  near  the  closed 
end  of  the  crack. 

Figure  23  describes  the  propellant  surface  temperature  distribution 
at  various  times.  The  curve  for  0.1  ms  shows  that  the  propellant  sur- 
face near  the  crack  entrance  is  slightly  above  the  ambient  temperature 
of  298°K  due  to  convective  heating  of  the  propellant  surface.  At  0.5 
ms,  the  propellant  surface  temperature  has  a sharp  gradient  near  the 
crack  entrance  caused  by  strong  convective  heating.  The  surface  temper- 
ature increases  continuously  until  the  state  of  ignition  is  reached.  In 
the  model,  the  surface  temperature  is  kept  constant  at  920°K  after  the 
propellant  is  ignited.  Therefore,  the  curves  for  the  times  of  1 and  6 
ms  have  flat  portions  near  the  entrance  of  the  crack  where  the  propel- 
lant is  burning.  The  propellant  surface  temperature  near  the  closed-end 
or  tip  of  the  crack  remains  low  since  convective  heating  of  the  propel- 
lant surface  is  significantly  lower  in  this  region  than  the  crack  en- 
trance region. 

Figure  24  shows  the  calculated  growth  of  the  crack  gap  due  to 
combustion.  Since  there  is  no  ignition  at  0.1  ms,  the  initial  crack  gap 
profile  is  the  same  as  the  curve  labeled  t » 0.1  ms.  For  the  time 
period  shown,  most  of  the  contour  change  occurs  near  the  open  end  of  the 
crack.  The  amount  of  surface  regression  is  the  highest  at  the  crack 
opening  since  the  propellant  has  been  burning  for  a longer  period  of 
time;  the  erosive  burning  effect  also  contributes  to  the  total  regres- 
sion in  this  region. 


PROBLEM  AREAS  IN  THE  STUDY  OF  PROPELLANT  CRACK  COMBUSTION 


In  recent  years,  a significant  amount  of  work  has  been  conducted  in 
the  investigation  of  crack  combustion  in  the  USA  and  also  in  other 
countries  such  as  the  USSR,  Japan,  and  England.  However,  the  convective 
burning  in  solid  propellant  cracks  is  not  yet  fully  understood  and  many 
problems  have  not  been  completely  solved.  It  is  beneficial,  in  the 
authors'  opinion,  to  list  some  of  the  unsolved  problems  here  so  that  the 
technological  gaps  in  this  area  are  identified.  It  should  be  noted  that 
some  of  the  problems  listed  are  the  immediate  subject  areas  that  are 
being  Investigated  at  the  Pennsylvania  State  University  but  other  prob- 
lems may  take  years  before  adequate  answers  are  obtained.  Hopefully  the 
following  list  of  questions  may  stimulate  further  Interest  among  the 
researchers  in  finding  solutions  to  these  problems. 
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FIGURE  22.  Calculated  Gas  Temperature  Distribution 
for  Various  Tiiues  (Simulation  of  DNC  Test  No.  21). 


NWC  TP  6049 


FIGURE  23.  Calculated  Propellant  Surface  Temperature  Distribu- 
tion at  Various  Times  (Simulation  of  DNC  Test  No.  21). 


FIGURE  24.  Calculated  Growth  of  Crack  Gap  Due 
to  Combustion  (Simulation  of  DNC  Test  No.  21). 


51 


NWC  TP  6049 


1.  Besides  the  dimensionless  groups  discussed  in  this  report, 
what  other  governing  parameters  are  important  to  consider  for  the 
coupled  solid-mechanic-combustion  problem  of  convective  burning  in  solid 
propellant  cracks? 

2.  When  the  hydraulic  radius  of  a crack  is  less  than  the  flame 
thickness  of  the  propellant,  the  combustion  process  is  incomplete  in  the 
local  region.  How  should  this  phenomenon  be  considered? 

3.  Is  dynamic  burning  significant  in  the  combustion  processes 
inside  a propellant  crack? 

4.  Belyaev,  et  al.  has  reported  that  the  velocity  of  the  flame 
front  depends  on  the  pore  shape.  For  example,  if  a flat  pore  is  re- 
placed by  a cylindrical  pore  with  the  same  hydraulic  diameter,  the  flame 
spreading  rate  is  much  lower.  Belyaev  stated  that  "there  is  no  satis- 
factory explanation  for  this  phenomenon."  The  parameters  such  as  ratio 
of  burning-to-wetted  perimeter,  surface  roughness  of  the  manufactured 
cracks,  and  secondary  flow  effects  may  have  a direct  influence  in  this 
phenomenon. 

5.  The  surface  roughness  of  a fissure  or  crack  in  a propellant 
grain  can  differ  appreciably  from  that  of  a manufactured  crack.  The 
projections  of  AP  crystals  or  binder  in  a rough  crack  would  be  easier  to 
ignite  and  would  result  in  a faster  flame  spreading  rate.  How  should 
one  manufacture  cracks  that  are  closer  to  the  real  crack  surface? 

6.  Do  we  have  sufficient  information  on  the  heat  transfer  and 
friction  correlations  which  are  required  physical  inputs  in  the  theoret- 
ical modeling  of  crack  combustion? 

7.  What  is  the  best  igniton  criterion  to  use? 

a.  Is  a surface  temperature  criterion  appropriate?  How  does 
the  critical  surface  temperature  vary  with  presure,  flow  velocity, 
surface  roughness,  presence  of  stress-induced  subsurface  vacuoles,  etc.? 

b.  Is  the  "Dipstick"  (subsurface  temperature)  ignition 
criterion  suitable? 

8.  A melt  or  liquid  layer  may  be  formed  at  the  propellant  surface 
during  the  combustion  of  HMX  (at  certain  pressures)  or  some  binders  such 
as  polyurethane.  The  melt  layer  may  cover  the  pore  to  reduce  the  pore 
size.  This  efiect  may  be  beneficial  in  retarding  the  transition  from 
normal  deflagration  to  convective  burning.  Can  this  effect  be  adequate- 
ly described  in  the  theoretical  modeling? 
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9.  According  to  Boggs,  et  al. , cracks  can  develop  during  the 
burning  of  HMX  crystals  due  to  steep  subsurface  thermal  gradients. 
Unlike  the  melt  layer  effect  described  above,  these  micropores  should 
contribute  to  the  transition  from  normal  deflagration  to  convective 
burning.  How  should  these  micropores  be  treated  in  the  theoretical 
modeling? 


10.  For  metallized  propellants,  what  is  the  role  of  the  metal 
particles  in  the  ignition  and  combustion  processes  inside  cracks — 
especially,  when  the  melted  particle  diameter  is  of  the  same  order  as 
the  hydraulic  diameter  of  the  crack? 

11.  Closure  of  a burning  propellant  crack  due  to  stress  wave 
interactions  may  cause  extinction.  What  degree  of  closure  is  required 
to  cause  extinction? 

12.  The  accelerated  convective  burning  inside  propellant  cracks 
can  cause  the  formation  of  shock-like  fronts  in  the  gas  phase.  What  is 
the  mechanism  of  transition  from  these  pressure  fronts  to  a shock  or 
detonation  wave  in  porous  propellants?  What  role  do  stress-induced 
vacuoles  play  in  the  transition  from  convective  burning  to  low  velocity 
detonation  (LVD)  or  detonation  (adiabatic  compression,  etc.)? 

13.  What  is  the  functional  form  of  the  propellant  stress-strain 
equation?  Is  a nonlinear  viscoelastic  model,  o - o(t,c,e,T,P,  etc.), 
absolutely  essential?  Are  there  sufficient  data  available  to  charac- 
terize the  propellant  constitutive  equation? 

14.  In  the  thermal  penetration  layer  of  burning  propellants,  are 
thermal  stresses  significant?  More  importantly,  is  the  propellant 
thermorheologically  simple  for  temperatures  approaching  800-1000 °K? 
Answers  to  these  questions  are  necessary  in  selecting  the  number  and 
size  of  finite  elements  in  the  stress  analysis. 


SUMMARY  OF  PROGRESS  AND  CONCLUSIONS 


1.  The  governing  equations  for  describing  the  heat  transfer  and 
gas  dynamic  processes  in  solid  propellant  cracks  have  been  nondimension- 
alized.  A set  of  governing  dimensionless  parameters  has  been  obtained 
to  help  future  data  correlation  and  analysis. 

2.  A number  of  modifications  have  been  made  in  the  crack  combus- 
tion program.  The  major  ones  are  given  in  the  following: 
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a.  In  order  to  obtain  more  information  about  the  in-depth 
thermal  profile  of  the  propellant,  the  one-dimensional  transient  heat- 
conduction  equation  was  solved  by  an  explicit,  variable-mesh  finite- 
difference  method.  The  propellant  surface  temperature  can  be  determined 
either  by  the  finite  difference  solution  of  the  partial  differential 
equation  for  heat  conduction  or  by  the  Goodman  integral  method. 7 The 
"Dipstick"  ignition  criterion,  postulated  by  Price*  has  been  incorpor- 
ated in  the  crack  combustion  program. 

b.  The  boundary  condition  at  the  entrance  of  the  crack  has 
been  improved.  The  program  can  now  more  accurately  describe  the  in- 
fluence of  rapidly  varying  rocket  chamber  pressurizations  on  the  subsonic 
and  sonic  inflow  of  gases. 

c.  The  restart  feature  of  the  crack  combustion  code  has  been 
improved  to  stop  the  calculation  at  any  desired  time  to  examine  the 
solution,  thus,  saving  both  time  and  money.  (The  most  up-to-date  version 
of  the  crack  combustion  program  was  sent  to  NWC  irt  February  1978.) 

3.  To  facilitate  the  photographic  recording  of  the  transient 
flame  spreading  event,  a reproducible  and  reliable  solid  propellant 
igniter  system  has  been  designed,  fabricated,  and  tested  to  replace  the 
previously  developed  gaseous  igniter.  The  purpose  of  this  replacement 
is  to  reduce  the  time  interval  between  the  onset  of  the  igniter  and  the 
first  ignition  inside  the  crack. 

4.  To  simulate  chamber  pressurization  during  the  ignition  trans- 
ient of  a rocket  motor  or  the  sudden  pressurization  for  an  initially 
submerged  crack  in  a solid  propellant,  the  pressurization  rate  of  the 
crack  combustion  chamber  is  controlled  by  varying  the  igniter  charge 
weight,  the  flow  areas  of  the  inlet  and  exit  nozzles  in  the  combsutlon 
chamber. 


5.  Numerous  test  firings  have  been  conducted  in  the  crack  combus- 
tion chamber  to  determine  the  influence  of  the  solid  propellant  igniter 
on  the  flame  spread  and  pressurization  in  single-pore  propellant  samples. 
The  previously  reported  debond  problem3  at  the  interface  between  the 
propellant  sample  and  the  brass  mold  has  been  solved. 

6.  A procedure  has  recently  been  developed  to  fabricate  pro- 
pellant samples  from  propellant  slabs  to  form  slit-type  cracks.  This 
in-house  fabrication  procedure  has  the  following  advantages : 

(1)  economy,  (2)  well  controlled  crack  geometry,  (3)  wide  range'  of  crack 
length  to  gap  width  ratios,  and  (4)  the  flexibility  of  viewing  the  flame 
spreading  from  different  orientations. 


Private  communication  with  Channon  F.  Price,  Naval  Weapons  Center, 
April  1977. 
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7.  The  crack  combustion  chamber  has  been  modified  and  various 
components  have  been  fabricated  to  accommodate  both  the  NWC  cast  speci- 
mens and  the  PSU  fabricated  propellant  cracks. 

8.  A literature  survey  on  crack  fracture  and  deformation  has  been 
conducted.  From  this  review,  a modified  version  of  NONSAP  (NFAP)  has 
been  selected  to  couple  with  the  crack  combustion  code  and  the  NFAP  code 
has  been  obtained  and  loaded  on  the  Pennsylvania  State  University  com- 
puter. 


9.  The  comparison  of  the  theoretical  results  with  experimental 
pressure-time  traces  showed  good  qualitative  agreement.  Although  the 
quantitative  comparison  showed  a reasonable  match,  some  differences  were 
noted.  In  order  to  improve  the  predictive  capabilities  of  the  theoreti- 
cal model,  the  following  physical  phenomena  must  be  considered: 

a.  incorporation  of  solid  mechanics  of  the  propellant  into 
the  crack  combustion  model; 

b.  possible  ignition  of  propellant  in  the  crack  tip  region 
caused  by  a mechanism  other  than  convective  heating  of 
the  propellant  surface;  and 

c.  dynamic  burning  behavior  of  the  propellant  under  rapid 
transient  pressure  excursions. 
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Appendix  A 

CHECK  LIST  OF  CRACK  COMBUSTION  TEST  PROCEDURE 


1.  Focus  the  camera  and  set  the  f-stop  at  the  desired  position. 

2.  Check  that  the  START/STOP  switch  on  the  camera  speed  controller  is 
in  the  STOP  position. 

3.  Make  sure  that  test  cell  power,  ignition,  and  sclenoid  de- 
activation switches  are  all  in  OFF  (down)  position  prior  to  film 
loading  in  camera. 

4.  Record  the  framing  rate  (PPS)  from  the  setting  on  speed  controller 
and  the  f-stop  from  camera. 

5.  Set  the  LED  driver  switch  to  the  ON  position.  Set  the  pulse  rate 
(pulses/s)  of  the  LED  maker  and  record  the  rate. 

6.  Check  to  ensure  that  the  camera  and  LED  driver  electrical  plugs  are 
securely  inserted  in  power  receptacles. 

7.  Mount  the  igniter  system  onto  the  crack  combustion  chamber.  Leave 
the  safety  plug  in  the  solenoid  armature. 

8.  Establish  flow  of  water  to  transducers  and  check  for  leaks. 

9.  Check  firing  relay  box  to  ensure  that 

a.  Common-time  leads  are  connected  to  both  the  camera  and  tape 
recorder  inputs 

b.  110  VAC  solenoid  lead  la  connected  to  relay  from  power  source 

c.  12  VDC  connected  to  electromagnet  Input 

10.  Activate  ONLY  test  cell  power  switch.  (Must  be  on  for  a minimum  of 
20  seconds  prior  to  firing  with  the  camera  speed  controller.) 

11.  Perform  a final  check  of  switches  and  leads 

a.  Transducers  leads,  charge  amplifiers,  and  tape  Inputs  conform 
to  that  previously  checked  on  DATA  SHEET 

b.  Charge  amplifiers  in:  (1)  CHARGE  mode,  (2)  GROUNDED  position 
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c.  Both  toggle  switches  on  calibration  panel  are  In  OFF  (up) 
position  and  rotary  switch  is  in  a null  position. 

12.  Turn  ON  test  cell  exhaust  fan. 

13.  Perforin  the  following  actions  in  the  test  cell: 

a.  Check  that  the  LED  driver  power  light  is  ON 

b.  Ensure  that  the  common-time  trip  indicator  on  LED  driver  is 
not  lit 

c.  Connect  lead  from  relay  box  to  solenoid  armature 

d.  Turn  ON  camera  lights 

e.  ARM  THE  SYSTEM  by  removing  the  solenoid  safety  plug 

f.  Close  the  test  cell  door 

14.  Firing  Officer  then: 

a.  Starts  tape  drive  in  RECORD  mode 

b.  Says:  "DNC  Test  No.  the  date,  "Unground  Charge  Ampli- 

fiers," the  "Fire"  (while  simultaneously  activating  the  IGNI- 
TION switch) 

c.  Immediately  after  percussion  element  discharges,  turn  ON  (flip 
up)  the  solenoid  de-activation  switch 

d.  NEVER  turn  off  test  cell  power  or  ignition  switches  while  the 
combustion  event  is  occurring  . Allow  ample  time  for  the 
entire  film  to  transport  to  take-up  reel. 

e.  STOP  tape  immediately  after  the  firing  is  completed  and  then 
GROUND  charge  amplifier. 

f.  Turn  OFF:  camera  lights;  test  cell  power,  ignition,  and 
solenoid  de-activation  switches;  and  then  cooling  water  to 
transducers  if  any  leaks  appear  in  water  lines. 

15.  Put  calibration  signals  on  magnetic  tape. 

a.  Start  tape  drive  in  RECORD  mode  and  say  "Calibration" 

b.  Simulate  a pressure  in  the  charge  amplifiers,  by  applying  a 
negative  DC  voltage  to  the  charge  amplifier  calibration  in- 
puts, then  switch  charge  amplifiers  to  both  CAL  mode  and 
OPERATE  position  and  then  switch  second  toggle  switch  on 
calibration  panel  down  for  say  only  3-5  seconds.  Finally, 

GROUND  charge  amplifiers  and  turn  off  the  tape  recorder. 

16.  Use  a black  hood  when  removing  film  from  the  camera  and  specify  all 
information  on  the  film  box. 

17.  Put  all  remarks  on  DATA  SHEET  however  trivial  they  may  be  (i.e., 
did  all  the  crack  burn,  etc.).  Use  a Kodak  Instamatic  or  press 
camera  to  record  all  interesting  disassembly  views. 
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Burning  surface  area 

Andreev  number  2 

Cross-sectional  area  of  the  crack,  cm 

The  spatial  change  of  cross-sectional  area  of  crack  with  respect 
to  the  axial  distance,  cm 
Area  of  inert  wall 

Pre-exponential  factor  in  the  nonerosive  burning  rate  law, 
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Hydraulic  diameter  of  the  crack,  cm 

Total  stored  energy  (internal  and  kinetic),  cal/g  2 
Acceleration  of  gravity,  conversion  factor,  g-cm/gf-sz 
Depth  of  the  gap  in  flow  visualization  tests,  cm 

2 

Local  convective  heat  transfer  coefficient,  cal/cm  -s-°K 

Local  convective  heat  transfer  coefficient  over  the  propellant 
surface,  cal/cm2-s-°K 

Local  convective  heat  transfer  coefficient  over  non-propellant 
pore  wall,  cal/cm2-s-°K 

Enthalpy  of  combustion  gas  at  adiabatic  flame  temperature,  cal/g 

Mechanical  equivalent  of  heat,  gf-cm/cal 

Erosive  burning  constant,  cm3-°K/cal 

Length  of  the  crack,  cm 

Molecular  weight,  gf/g-mole 

Pressure  exponent  in  the  nonerosive  burning  rate  law 

Static  pressure,  g^/cm2 
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Burning  perimeter,  cm 

wetted  perimeter  of  the  port,  cm 

Specific  gas  constant  for  the  combustion  gases,  g^-cm/g-°K 
Reynolds  number 

Burning  rate  of  the  solid  propellant,  including  the  erosive 
burning  contribution,  cm/s 
Stanton  number 

Temperature  (without  subscript,  static  gas  temperature),  °K 

Average  film  gas  temperature  (T  + T )/2,  °K 
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Tf  Adiabatic  flame  temperature  of  the  solid  propellant,  °K 

Tpi  Initial  propellant  temperature,  °K 

Tps  Propellant  surface  temperature,  °K 

1'ws  Nonpropellant  wall  surface  temperature,  °K 

t Time , s 

u Gas  velocity,  cm/ s 

Vgf  The  velocity  of  propellant  gas  at  the  burning  surface,  cm/s 
The  coefficient  in  viscosity-temperature  relation 
x Axial  distance  from  propellant  crack  opening,  cm 
xl  Position  at  the  end  of  crack,  cm 

Xp  Axial  distance  along  the  crack  at  which  propellant  begins,  cm 
y Perpendicular  distance  from  the  propellant  surface  into  the 
solid,  cm  2 

a Thermal  diffusivity,  cm  /s 
3 Erosive  burning  exponent 
Y Ratio  of  specific  heats  * 

6 Gap  width  of  the  crack,  cm 
e A small  number  of  °K 

e Surface  roughness,  cm 

0 Angle  measured,  in  a counterclockwise  direction,  at  the  lower  side 
of  the  propellant,  degree 
X Thermal  conductivity,  cal/cm-s-°K 

p Gas  viscosity,  g/cm-s  (poise) 

p Density  (without  subscript,  gas  density),  g/cm 
Shear  stress  on  the  port  wall,  g^/cm^ 

Normal  viscous  stress,  gf/cm^ 

Dimensionless  parameters 
* Reference  parameters 


Subscripts 


c Rocket  chamber 

cri  Critical  condition  for  surface  ablation 

fp  Flame  propagation 

fs  Flame  stand-off 

i Initial  value 

ign  Ignition  condition 

pr  Propellant 

th  Thermal  wave  thickness 


:emperature) , °K 
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